Enhancement of charged macromolecule capture by nanopores in a salt gradient 
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Nanopores spanning synthetic membranes have been used as key components in proof-of-principle 
nanofluidic applications, particularly those involving manipulation of biomolecules or sequencing of 
DNA. The only practical way of manipulating charged macromolecules near nanopores is through 
a voltage difference applied across the nanopore-spanning membrane. However, recent experiments 
have shown that salt concentration gradients applied across nanopores can also dramatically en- 
hance charged particle capture from a low concentration reservoir of charged molecules at one end 
of the nanopore. This puzzling effect has hitherto eluded a physically consistent theoretical expla- 
nation. Here, we propose an electrokinetic mechanism of this enhanced capture that relies on the 
electrostatic potential near the pore mouth. For long pores with diameter much greater than the 
local screening length, we obtain accurate analytic expressions showing how salt gradients control 
the local conductivity which can lead to increased local electrostatic potentials and charged analyte 
capture rates. We also find that the attractive electrostatic potential may be balanced by an out- 
ward, repulsive electroosmotic flow (EOF) that can in certain cases conspire with the salt gradient 
to further enhance the analyte capture rate. 

PACS numbers: 



INTRODUCTION 

Recent interest in electrokinetic manipulation of 
charged macromolecules has been motivated by techno- 
logical applications, particularly those involving sorting 
and sequencing of nucleic acids. In a typical realization of 
single molecule DNA sequencing, an ionic solution is sep- 
arated by a membrane with a small pore across the mem- 
brane, connecting two otherwise separated bulk reservoirs 
(cf. Fig. 1). When an electric potential is applied across 
the membrane, ionic current flowing through the pore is 
detected. DNA and protein molecules placed on one side 
of the membrane (the right reservoir in Fig. 1), even 
at low concentrations, can occasionally block the pore, 
reducing the ionic current. A time trace of the ionic cur- 
rent flowing across the pore therefore directly measures 
the statistics of blocking and unblocking events. 

Experimentally, numerous modifications of the basic 
configuration have been studied. Biological pores such 
as a-hemolysin have also been chemically modified to al- 
ter internal charges, leading to possible enhancements of 
the capture frequency and translocation rates of biopoly- 
mers through the pore [1]. A number of groups have 
also recently fabricated synthetic pores [2-4], typically 
through SiN membranes, for use in macromolecule cap- 
ture experiments. Besides pore design, other approaches 
to better control macromolecular analyte (both charged 
and uncharged) capture and translocation have been ex- 
plored. In recent measurements using synthetic pores, 
an enhanced capture rate of DNA was observed in the 
presence of a salt gradient [5, 6]. Salt added to the op- 
posing, non-analytc reservoir increased the capture rate 
as a function of the salt ratio between the two reservoirs. 

Much of the theoretical effort, including molecular sim- 



ulations, has concentrated on the physics of polymer 
translocation through the pore [7-9]. However, since 
macromolecular capture is sensitive to applied voltage, 
occurs even with uncharged molecules, and is a stochas- 
tic process, the relevant mechanism will involve the in- 
terplay among electrostatics, fluid flow, and the statistics 
of capture. Although the electric potentials and flows in- 
side a nanopore have been presented in the context of 
Poisson-Nernst-Planck models [10, 11], the dynamics of 
initial particle capture requires a more detailed analy- 
sis of the field configuration in the bulk reservoir, near 
the pore opening. Finally, the analyte density in the 
reservoir needs to be determined [2] in order to solve the 
problem of capture of a charged macromolecule to the 
pore mouth in the presence of fluid flow and electroki- 
netic forces that stochastically switch according to pore 
blockage. The polymer capture problem has been treated 
by Wong and Muthukumar [12], who derive capture rates 
in the presence of electroosmotic flow (EOF) induced by 
an applied voltage bias and surface charges on the inner 
pore surface. In this study, as in many others, [8, 12, 13] 
the electric field in the reservoir was neglected and the 
capture problem was solved only in the fast translocation 
limit. 

In this paper, we show that bulk electric fields and 
a more detailed calculation of the capture problem are 
necessary to explain the recently-observed salt gradient- 
induced enhancement of charged analyte capture rates. 
By using geometrical approximations in the high salt 
limit, we find analytic expressions for the local salt con- 
centration and electrostatic potential that are accurate 
for a range of parameters relevant to typical experiments. 
In the presence of an imposed salt concentration gradi- 
ent, we show that although most of the potential drop 
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occurs across thin pores, the electrostatic potential at 
the pore mouth plays an important role in the analyte 
capture rate and cannot be neglected. We also show how 
an EOF into the analyte reservoir (by itself convecting 
macromolecules away from the pore) can interact nonlin- 
early with the salt gradient and electrostatic potential to 
actually increase the capture rate. Finally, we analyze 
implicit solutions for mean capture rates using the full 
steady-state Debye-Smulochowski capture problem with 
a partially absorbing boundary for the analyte concentra- 
tion at the pore mouth. 




right reservoir 



charged 
macromolecular 
analyte 




ELECTROKINETIC EQUATIONS 

A typical experimental set-up is depicted in Fig. 1(a) 
where two reservoirs containing aqueous solution are sep- 
arated by an electrically insulating membrane contain- 
ing a single conducting nanopore through which fluid can 
flow. A voltage bias v is applied across electrodes placed 
far from the pore, resulting in an ionic current through it. 
The full steady-state electrokinetic equations for the local 
electrostatic potential <j>(r,t), ion concentration d(r,t), 
and local fluid velocity U(r, f): 



V • (eVp(r,*)) = -47re^z lCl (r) 



(1) 



u(r) • Vc 4 = V • (DiVci) + e^-^-V • (aV<p). (2) 



and 



= -Vp + r/V 2 u + eVipY^ z lCl {r), 



(3) 



where Zi is the valency of solute species i, e is the dielectric 
permittivity of the solution, Di is the diffusivity of species 
i, i] is the dynamic viscosity of the solution, and p(r) is 
the local hydrostatic pressure in the fluid. 

For simplicity, we henceforth consider a two component 
ionic solution with z± = ±1 where both ion species have 
the same diffusivity D + = D- = D lon . Reactions at the 
electrodes will also be symmetric such that no net charge 
is built up in the bulk solution. We analyze Eqs. 1-3 
in the geometry shown in Fig. 1 where, for mathemati- 
cal convenience, the electrodes far from the nanopore are 
assumed to be hemispherical caps with radius r — ► oo. 

The boundary conditions far away from the nanopore 
at these "infinite" electrodes arc c + = c_ = Cj,, u = 0, 
and if = v in the left reservoir, and c+ = c_ = cr, u = 0, 
and if = in the right reservoir. To obtain effective 
boundary conditions at the membrane and on the inner 
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FIG. 1: Schematic of electrokinetic focussing experiments, (a) 
Charged analytes are placed in the right reservoir and a volt- 
age bias is applied, (b) The structure of the pore, electrostatic 
potential, and electroosmotically driven fluid flow in dimen- 
sionless units, with distance measured in units of the pore 
radius. For long pores (a/d = w <C 1), the concentration 
and potential fields are approximated as constant within the 
small hemispherical regions capping the pore (denoted by ±). 
Field and flux continuity conditions are applied at the hemi- 
spherical surfaces. The electroosmotic flow field in the right 
chamber will be approximately spherically symmetric (red) if 
the membrane surface is uniformly charged but will be more 
lobe-like (black) if the membrane flange is uncharged and no- 
slip boundary conditions are imposed [14, 15]. 



surface of the pore, we first extract the appropriate physi- 
cal limit by defining distance in units of the pore diameter 
a, and nondimcnsionalizing Eqs. 1-3 according to 



r - c± 

CR 



$ = 



etp 



v = 



u = 



-u, P = 



P 



An, 



k B Tc R 



M = 



'/ 



(4) 



ksTcRa 2 



Equation 1, the sum and difference of the z = +1 and z = 
— 1 components of Eq. 2, and Eq. 3 become, respectively, 



V-( £ V$(r))+A J? Q(r)=0 (5) 

V • [VE(r) + Q(r)V$(r) - U(r)E(r)] = (6) 

V • [VQ(r) + E(r)V$(r) - U(r)Q(r)] = (7) 

AiV 2 U(r) - VP(r) +Q(r)V$(r) = 0. (8) 

Above, we define the ion concentration difference Q(r) = 
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±(C+(r) - C_(r)), sum E(r) = i(C+(r) +C_(r)), and 
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87re 2 c^a 2 
k B T 



{k R o) 2 . 



(9) 



The quantity represents the inverse ionic screening 
length associated with the ionic solution deep in the right 
reservoir. For sufficiently high salt concentrations cr ~ 
1 — 0.1M, the corresponding screening length k^ 1 ~ 0.3 — 
lnm, while synthetic pores have radii on order of at least 
a = 3 — 5mn, rendering A > 10 large. At distances of 
least a screening length away from the membrane or pore 
surfaces, this separation of scales allows us to consider 
only the "outer" solutions associated with charge-neutral 
conducting fluid bulk reservoirs where boundary layers of 
charge separation do not reach. 

In the A^ 1 = e — > limit, Eq. 5 represents a singu- 
lar perturbation. The outer solutions to the system of 
equations can be found by considering expansions of the 
form 



(10) 



Upon using these expansions in Poisson's equation, eV • 
(eV$) = Q, we find an outer solution Qo = 0, as expected 
in the charge-neutral limit of a fluid conductor. To find 
solutions accurate to 0(e ), we must solve the remaining 
equations 
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P(r) 
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[VE (r)-U (r)E (r)] = 
V ■ [E (r)V$ (r)] = 
M V 2 U (r) - VP (r) = 



(11) 
(12) 
(13) 



Henceforth, we consider only the zeroth order solutions 
and drop the subscript (0) notation. Equation 11 simply 
describes the steady-state convection-diffusion of the to- 
tal local salt concentration E(r). The ionic conductivity 
of a locally neutral electrolyte solution is proportional to 
this local salt concentration. Eq. 12 is an expression of 
Ohm's law that includes a spatially varying ionic conduc- 
tivity. 

In the relevant limit of small screening length, the near- 
field boundary conditions for the concentration and po- 
tentials associated with the outer Eqs. 11-12 can be ob- 
tained by noting that the inner solution near the mem- 
brane decays exponentially {e.g., as does the solution 
to the Poisson-Boltzmann equation). To match these 
exponentially decaying inner solutions requires effective 
Neumann boundary conditions for the outer solutions of 
E and <J> near the solution-membrane interface. Such 



boundary conditions embody the impenetrable nature of 
the membrane to both the salt and the ionic current, and 
has been more formally derived in [10]. The Neumann 
boundary conditions allow the use of spherical symme- 
try of the solutions within each reservoir, rendering Eqs. 
11 and 12 one-dimensional in the radial coordinate. The 
far-field boundary conditions on the dimensionless quan- 
tities are E r (r — > oo) = 1, E^(r — > oo) = E/,, and 
$c(r — * oo) = V, & r {r — ► oo) = 0. Although complicated 
expressions have been derived to compute the exact con- 
centration profiles associated with steady-state diffusion 
through a finite width pore [16], our spherical approxi- 
mation dramatically simplifies the analysis of Eqs. 11-12. 

Finally, we must also consider the structure of the flow 
field U in Eq. 13, even if no hydrostatic pressure gradient 
is externally applied and P{r — ► oo) = 0. A nonzero outer 
velocity field will arise from boundary charge-induced 
electroosmosis, in which a surface charge on the inner sur- 
face of the pore enhances the concentration of the countc- 
rions near the inner surface the pore. An applied electric 
then pushes this slightly charged layer, dragging the bulk 
fluid with it. Matching the outer velocity field with the 
electric field-driven, inner layer of fluid within a Dcbye 
screening length of the charged inner pore surface results 
in a plug-like outer flow velocity profile. The plug-like 
profile allows us to consider the outer solution for U in- 
side the pore as a constant Ux. The magnitude U of 
the plug flow velocity is proportional to the "^-potential" 
at the pore surface and the electric field applied across 
the pore [17]. Within linearized electrostatics described 
by the Debye-Huckel equation, the ^-potential is propor- 
tional to the surface charge density and the local screen- 
ing length of the solution. 

In the absence of pondcrmotive body forces on the 
fluid outside the pore, the velocity field U, under no- 
slip boundary conditions, has been calculated as a series 
expansion [14] and in terms of integrals over Bessel's func- 
tions [15]. The lobe-like flow patterns (cf. Fig. 1(b)) were 
obtained using no-slip boundary conditions at the mem- 
brane interface and are not spherically symmetric. How- 
ever, if the membrane walls are also uniformly charged, 
potential gradients along the wall exterior to the pore 
will drive an inner EOF along the wall, giving rise to an 
effective slip boundary condition for the outer flow field 
U near the wall. We will show that the outer-solution 
electric field tangential along the membrane falls off as 
1/r 2 , plus logarithmic corrections in the presence of salt 
gradients. Therefore, the velocity profile in the reservoirs 
can be self-consistently approximated by a radial profile 
U = Uf/r 2 obeying incompressibility. The actual veloc- 
ity profile will qualitatively resemble radial flow, while 
retaining some features of a lobed-flow profile. Since we 
expect our gross results will be sensitive mostly to the 
typical magnitude of U, and not to the details of its an- 
gular dependence, we will also assume a simple radial 
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form for U(r). 

Since we assume radial symmetry in the outer solutions 
of all quantities, we apply uniform boundary conditions 
at both the near and far boundaries at r — ► oo and r = 1, 
respectively. Upon defining the variables 



G = cxp 



and H : 





~Ud~ 




~U~ 


exp 




= cxp 






a 




w 



(14) 



where w = a/d <C 1 is the pore aspect ratio, the total ion 
concentration can be found in terms of the normalized 
bulk salt concentration E^ far from the pore in the left 
reservoir. When the pore is unblocked by analyte, salt 
freely diffuses across. Using the general solutions of Eq. 
11 in each region (left and right reservoirs, and pore), 
and imposing conservation of ion flux at the surfaces of 
the hemispherical caps between these regions (cf. Fig. 
1(b)), 2ird r E e (r = 1) + nd x E p (x = -l/(2w)) = and 
27r<9 r .£ r (r = 1) - nd x H p (x = l/(2w)) = 0, we find 



Mr) = 



H L GH-1 (Si-1) 



GH — 1 



GH-l 



exp 



U 
27 



(15) 



in 



= ($+-$_)- 



1 + 2w 



hi 



<*>-, 



and 



/r(y a E r (y))-id|, 

ln£ r (r) + U/{2r) 
lnE+ + U/2 ■ 



(19) 



(20) 



In order to explicitly determine the potentials <I>±, wc 
apply Kirchhoffs Law, 

27rI]_cV$<> (r = 1) + nE-dx^pix = -d/(2a)) = 0, 
27rE + a r $ r (r = 1) - 7rZ + d x <S> p (x = d/(2a)) = 

(21) 

conserving ionic current across the pore mouths to find 



T, p (x) 



E L GH-1 {E L -1)VGH 



GH-l 



GH-l 



(16) 



$+ =V 



ln£+ + U/2 
lnE L + £/(l + l/w) 



(22) 



and 



and 



Er(r) 



Y^lGH-I 
GH - 1 



1)GH 



GH-l 



■ cxp 



U_ 

2r 



(17) 



Although exact solutions can be computed [16], for small 
the simple forms given in Eqs. (15-17) are accu- 
rate to 0(w). 

The salt concentration £(r) determines the local ionic 
conductivity as a function of the EOF velocity U and the 
experimentally imposed bulk salt ratio E^. We now sub- 
stitute the approximate solutions for E(r) into Eq. 12 
to find the electrostatic potential <f>(r). Denoting the ion 
concentration and electrostatic potential within the hemi- 
spheres capping the pore as E p (l/(2w)) = E r (l) = £ + , 
E p (l/(2w)) = E<(1) = and %{l/{2w)) = $ r (l) = 
$ p (— l/(2u;)) = $*(1) = and applying all far- 
field boundary conditions, we find 



$_ = V 



ln£_ +U(l + 2/w)/2 



lnE L + [/(l + l/w) 
In the limit of vanishing salt gradient, E^ — > 1, and 



$+(E L -> l,w,t/) 



wU 



(23) 



0(£ L - 1), (24) 



2(1 + to) 

while for vanishing clcctroosmotic flow, C/ — s- and 



$+(E L ,™,[/^0) 



V 



InE, 



In 



2 + ui(E L + 1) 



2(w + 1) 



(25) 



In the case E^ — > oo, the negligible resistance offered 
by the left reservoir makes <&+ slowly (logarithmically) 
approach V. 



Mr) =(V 



In 



V- 



2 



J7 
2r 



hi 



lnffH+f 



ln(|^ 



J7 

2r 



(18) 



RESULTS AND ANALYSIS 
Concentrations and potentials 

Fig. 2(a) shows the total salt concentration along the 
axial coordinate, made up of solutions to E^(r), T, p (x), 
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Effect of electroosmotic flows 

Now consider how the details of the electroosmotic flow 
velocity U may affect the electrostatic potential. For an 
EOF with given magnitude U, the potentials <I>± can be 
calculated from Eqs. 22 and 23. However, recall that the 
EOF velocity is driven by the potential difference applied 
across the pore [17]; therefore, U must be self-consistently 
solved by finding the root to 



FIG. 2: Salt concentration and electrostatic potential across 
an unimpeded pore of aspect ratio w = a/d — 0.1. (a) The salt 
density E(£) as a function of the axial coordinate for various 
EOF velocities U and salt ratio El = 3. (b) The normalized 
potential $(£)/V for various salt ratios El. 



and E r (r). For zero EOF (U — 0), the total salt concen- 
tration is linear in the pore region, but has a decaying 
structure in the two reservoirs. 

The spatially varying total salt concentration results in 
a spatially varying conductivity. Fig. 2(b) shows the re- 
sulting electrostatic potential along the axial coordinate 
£. For vanishing salt gradient El = 1, Eq. 20 shows the 
potential decays as l/r in the analyte reservoir, consis- 
tent with our assumption of a hemispherical EOF profile 
far from the pore. In the presence of an imposed salt 
gradient (El ^ 1), additional logarithmic terms arise in 
the reservoir potential $ r (r). These results show that 
although most of the potential drop occurs across the 
membrane-spanning pore, in the presence of ionic cur- 
rent, the electrostatic potentials in the reservoirs vary 
slowly with distance from the pore mouths. 

For Ex, = 1 and long pores (w <C 1), the poten- 
tial at the pore mouth in the grounded analyte reservoir 
~ wV may not be small for sufficiently large applied 
bias voltage. Furthermore, we will show that the analyte 
capture rate is sensitive to $ + as it depends on the ex- 
ponent of $ + . For imposed salt gradients El > 1, the 
conductivity in the left reservoir is relatively higher and 
more of the voltage drop occurs across the pore and the 
analyte (right) reservoir. This further raises the poten- 
tial $ r (r = 1) = $ + felt by the charged macromolecules 
at the mouth of the pore in the analyte chamber. Note 
that the potential $ r (r-) depends on the gradient of con- 
ductivity (arising from the salt concentration gradient), 
and not on its absolute value. The conductance vari- 
ation arising from salt concentration gradients provides 
a simple mechanism by which the potential $ r (r) can 
be increased through the salt ratio El, enhancing in the 
capture rate of charged analytes. 



U = T(U, E L )(p+[U, E L ] - *-[U, E L ]), 



(26) 



where we have explicitly denoted the dependence of <I>± 
on U and E^. The prefactor T(U, El) measures the 
effective electroosmotic permeability, which is inversely 
proportional to the pore length and fluid viscosity, and 
proportional to the "C-potcntial" [17]. Within the lin- 
earized Debye-Hiickel theory for electrolytes, this local 
^-potential is proportional to the pore surface charge 
times the local screening length k _1 (.t). In our prob- 
lem where the ionic strength is varying in the axial di- 
rection along the pore, the potential is also varying 
along the pore. Although a nonuniform surface poten- 
tial, along with the constraint of fluid incompressibility, 
can give rise to nonuniform flow within the pore, it has 
been shown that the net fluid flow across a pore can be 
found by averaging the ^-potential (or screening length) 
across the pore [18-20]: 



T(U,E L ) = 



ak B T\ 1 f d/2 dx 



ride J d 

• 1/(2 



Y R w 



d/2 K(x) 

d.r 



(27) 



-l/(2ti>) V%W 



where 



eac R 



A V2 



(28) 



is the dimcnsionlcss permeability referenced to the salt in 
the right chamber. For a typical SiN membrane with con- 
stant surface charge density \<j\ < 1 microcoulomb/cm 2 
[21], a 5nm radius, 40nm length pore in water gives 

|r«| < 0.2. 

Figure 3(a) shows the self-consistent EOF velocity U 
(obtained by using Eq. 27 and solving Eq. 26) as a func- 
tion of T/j, for various salt ratios El- For > and 
U > 0, salt is being swept from the left to right reservoirs. 
When El > 1, the pore feels a higher averaged salt con- 
centration, lowering the ^-potential, thereby reducing the 
EOF response to Tr. Conversely, when El < 1, a higher 
^-potential and stronger response arises. Fig. 3(b) shows 
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y=10, w=a/d = 0.2 IfelO, w=a/d=0.2 F=10, w=a/tf = 0.2 V=\0,w=a/d = 0.2 




FIG. 3: (a) EOF velocity U as a function of effective pore EOF 
permeability Fr. The response deviates from linear for large 
and small salt ratios. The deviations are most pronounced for 
Fr > where EOF is into the right reservoir (U > 0). (b) 
EOF velocity as a function of salt ratio for various effective 
pore surface charge densities. For Fr > 0, increasing the 
salt ratio decreases the effective screening length in the pore, 
reducing the EOF velocity (cf. Eq. 26). When F R < 0, 
the salt in the right reservoir is swept into the left reservoir, 
keeping the screening length approximately k^ 1 throughout 
the pore and very little dependence on El arises. In both 
plots V = 10 and w = a/d = 0.2. 



the EOF velocity U as a function of salt ratio E l for vari- 
ous r^. For El — > oo, U vanishes along with the effective 
pore ^-potential. 

The pore mouth potential <& + felt by the charged 
analyte is shown in Fig. 4. As a function of pore 
charge/permeability T^, the potential $ + exhibits a max- 
imum(minimum) for El > 1(El < 1). This nonmono- 
tonic dependence arises because for El > 1 and small 
U > 0, the EOF changes the conductance structure so 
that <I> + initially increases. In other words, the largest 
voltage drop across the analyte (right reservoir) reservoir 
occurs at small, positive U. For larger [/, high salt is 
swept well into the right reservoir reducing the effective 
relative conductivity across the pore. Most of the voltage 
drop then occurs in regions away from the pore well in the 
right reservoir, diminishing <&_)_ to the value wV/(2 + 2w) 
expected in the uniform salt (El = 1) limit. For El < 1, 
more voltage drop occurs across the left reservoir, reduc- 
ing $ + below wV/(2 + 2w). Fig. 4(b) shows that as the 
pore potential <J>+ increases with the salt ratio El. This 
is consistent with the observed increase in analyte cap- 
ture rate with increasing salt in the analyte- free chamber. 
Also consistent with Fig. 4(a) is the slight increase, then 
decrease in $ + as is increased 



FIG. 4: (a). The electrostatic potential in the right reservoir 
(normalized by the applied potential V), as a function of salt 
asymmetry El for various pore EOF permeabilities r_R. Here 
the pore aspect ratio was set to w = a/d = 0.2. (b) The 
magnitude of the potential near the mouth pore as a function 
of EOF permeability Fr for various salt ratios. 



nanopore as functions of parameters such as the applied 
salt ratio El, applied voltage bias V, and and pore sur- 
face charge/permeability Tr. When a charged analyte 
molecule of size of order a enters the hemispherical cap, 
it blocks the pore and prevents ion transport. Such a 
nonspccifically adsorbed particle can spontaneously des- 
orb from the mouth of the pore with rate k g. Alterna- 
tively, as in the case of DNA, it may translocate through 
the pore to the opposite, receiving reservoir. Although 
translocation of a polymer involves many stochastic de- 
grees of freedom, we will lump the process into a single, 
effective rate fct, such that k^ 1 represents the typical time 
for the macromolecule to fully tunnel across the pore, al- 
lowing ionic current to flow again. 

If the hemispherical cap can accommodate at most one 
blocking macromolecule, its mean field, steady-state oc- 
cupation < 9 < 1 is balanced according to: 



k on p(l)(l-6) = (fc off + k t ) 



(29) 



where p(l) = lim r _ >1 + p(r) is the analyte concentration 
just outside the cap region, and k on is its adsorption rate 
into the hemisphere. To explicitly determine 9 we need 
to relate the unmeasurable, kinetically-determined p(l) 
with the experimentally imposed macromolecular density 
Poo = p( r ~ * oo)- This relationship is obtained by solving 
the mean field, steady-state convection-diffusion equation 
for the density in the bulk region: 



Charged analyte distribution 

We now model how both the approximate EOF veloc- 
ity fU/(2r 2 ) and the electrostatic potentials (Eqs. 18, 
20, and 19) affect the capture of charged analytes to the 



(l-0)V-[A(r)p(r)]=V», r > 1, (30) 



where 



A(r) = A{r)f 



U 



2Dr 2 



' dr 



(31) 
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is the normalized drift arising from both a hydrodynamic 
flow and an electrostatic potential induced by electroos- 
mosis and ionic conduction, respectively. The drift in- 
duced by the electrostatic potential is proportional to the 
effective number of electron charges q of the analyte. The 
factor (1 — 0) reflects our assumption that the EOF flow 
and ionic current is completely shut off when a macro- 
molecule occupies the cap (8 = 1). In writing Eqs. 29 and 
30, we have implicitly assumed that the time r e ff required 
to form the effective potential U/(2Dr) + qd r $ r (r) is 
much less than the typical times associated with dissoci- 
ation, association, or translocation: r c ff <C k~g, k~^, k^. 
For relevant parameters, T c g ~ a 2 /(ArD) ~ 10~ 10 s. Pro- 
vided all the rates are slower than the rate of reestablish- 
ing the effective potential, we can assume that the field 



and flow configurations instantaneously follow those cor- 
responding to whether or not the pore is blocked. 

The boundary condition associated with Eq. 30 is ap- 
plied just outside the hemispherical cap (r = I + ) and is 
determined by macromolecular flux balance through the 
interface [22, 23] 



Dd rP (r)\ r=1+ - (1 - 6) [f + DqV<P + ] p(l) 

= konp(l)(l-0)-k o s6. 



(32) 



Upon solving Eq. 30 for p(r) and using Eq. 29 for p(l), 
we find 



J 



p(r) = e (1 " 9) X A (y) d v 



he 

D 



2T 2 e- (1 



-9) fy A(y')dy' dy . 



(fc off +kt)0 

Mi - 6 ) 



(33) 



r 



After setting r — > oo in Eq. 33, we relate the mean cap 
occupation 9 to the bulk analyte density p^ through the 
physical solution of the integro-transcendental equation 



W r/m , (k s + h)0 



Xl-B)(U/(2D)-q* + ) 



(34) 



where 



1(9) 



»(l-«)(5gr-9»r(r)) 



(35) 



When fc t = 0, adsorbed macromolecules do not translo- 
cate and can only detach back into the right bulk reser- 
voir. In this limit, 1(9) does not arise in Eq. 34 and 
8 depends only on the value <f> + of electrostatic po- 
tential at the pore mouth. Moreover, the dependence 
on the normalized macromolecule diffusivity arises only 
in U/(2D), the drift due to EOF at the pore mouth. 
Although the EOF and ionic conduction in the fluid 
arises from noncquilibrium processes, the density pro- 
file p(r) in the fct = limit is an equilibrium den- 
sity self-consistcntly determined by the effective poten- 
tial (1 — 9) [U/(2Dr) — q® r (r)]. For parameters such that 
A(l) ee [U/(2D) - Q$+] > ln( Poo fc on /fc off ), 



fcoff 



-C//(2D)+g* 4 



< 1. 



(36) 



This form corresponds to an equilibrium "adsorption 
isotherm" on the single site and depends only on the value 
of the potential energy at that site. 



Conversely, in the limit of high translocation rates 
such that k t k on /D > k off and [U / (2D) - q$ + ] > 
ln^oofcon/fct), the first iteration about 9 = to Eq. 34 
yields 



(k o Jk t )D Poo e- u /(2 D )+0*+ 
fc on /(0)e- c/ /( 2 - D )+9*+ + D 



(37) 



Here, the surface density p(l) ~ kt9/k on resembles that 
of an adsorbing sphere with an attachment rate fc on , but 
modified by the 0(1) term /(0) e - c/ /( 2I5 ' + « <E, + resulting 
from stochastically switching of the effective drift. 

Although the full solution for the pore occupation must 
be solved numerically, we see that 9 is exponentially sen- 
sitive to both the magnitude of the EOF velocity U and 
the electrostatic potential at the cap 4> + through an ef- 
fective drift defined by the combination U/(2D) — q&+. 
Note that both U and <f> + depend linearly on the bias 
voltage V, but nonlinearly on the salt ratio S^. The 
EOF velocity U is a function of through the solution 
Eq. 26, while depends indirectly on Tr through the 
resulting flow velocity U that changes the local conduc- 
tivity when Yil ^ 1. However, only the electrostatic drift 
g<3?+ depends on the effective analyte charge q. 

Since many analyte capture experiments exhibit infre- 
quent pore blocking, even with bias voltages as high as 
+250mV (V w +10), we use parameters that yield small 
numerical values of 9. Henceforth, we set the relative ana- 
lyte/ion diffusivity D = 0.01, p^ = 1Q~ 7 (corresponding 
to an analyte concentration of ~ InM for a ~ 5nm), and 
the effective analyte charge q — 30 (corresponding to that 
of an approximately 500bp strand of dsDNA [24] ) . Figure 
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V=W,w=a/d = 02 r„=0.1, w=aJd = 0.2 




FIG. 5: (a) Pore mouth occupation fraction 9 as a function of 
salt ratio El for various EOF permeability Fr and V — 10. 
Note the sharp increase in 9 as a function of El for positive 
surface charges. Parameters used were p x = 10~ 7 , V = 10, 
w = 0.2, k on /k ff — 1000, and k on /k t = 10 4 . For small El, 
larger Fu > induces larger U > 0, pushing the analyte away. 
At higher salt El, the EOF is mitigated due to the reduction 
in ("-potential (or effective surface charge) indicated by Eq. 27. 
The reduction of EOF to modest values allows the attraction 
from the term g<3>+ to overcome the repulsive effect of the 
EOF (cf. Fig. 4(a)), increasing 9. (b) Occupation fraction 
9 as a function of bias voltage V at different salt ratios and 
F R = 0.1. 

5 shows representative numerical solutions of Eq. 34 as a 
function of (a) salt ratio El, and (b) bias voltage V. Fig. 
5(a) shows that for larger r# (larger U), the repulsive 
U/(2D) term dominates in keeping 9 small. However, 
upon increasing El, attraction arising from an increas- 
ing term eventually increases 9. Larger values of 
El also attenuate the pore-averaged ("-potential, reduc- 
ing the repulsive EOF, further increasing 9. For fixed V, 
we will show that the analyte capture rate is proportional 
to 9; therefore, Fig. 5(a) predicts the capture rate as a 
function of salt ratio. 

To determine how 9 depends on the applied voltage, we 
use simple assumptions to approximate how the kinetic 
rates depend on V , 

fe off = L} ofi e~ fqV , k on = Won, h = w t V, (38) 

where u> g,u> on arc intrinsic detachment and attachment 
rates, and 1/uJt is the typical conditional mean time 
for analyte translocation across the nanopore under a 
v ps ksT voltage bias (translocation is assumed negli- 
gible when V = [1] but can be approximated for poly- 
mer translocation [25, 26]). When the macromolccular 
analyte blocks the pore and = 1, a fraction / of the 
q charges may be exposed to the potential in the pore. 
When the pore is completely blocked, this potential is 
approximately V since there is no voltage drop across 
the left reservoir and the pore. For detachment to oc- 
cur, an energy barrier fqV must be overcome, resulting 
in fe off ps uj oS e~ fqV [27]. In Fig. 5(b), 9 is plotted as a 



function of voltage V for various fixed salt ratios El- The 
kinetic parameters chosen were / = 0.02, ui s = uJ on , and 
Woff/wt = 10 5 . This choice of intrinsic rates corresponds 
to approximately 4% of captured analyte being translo- 
cated at V = 10, and 96% detaching back into the bulk 
analyte reservoir. For the chosen parameters, increasing 
a small voltage V raises 9 exponentially if there is an 
appreciable ratio El that enhances the positive EOF ve- 
locity U to increase beyond the repulsion caused by 
the positive flow U. At larger V, the larger EOF velocity 
not only pushes the analyte faster from the pore, but also 
contributes to reduction in the potential <f> + , resulting in 
a slower increase in occupation. 

Analyte capture rates 

We now define the mean analyte capture rates mea- 
sured in experiments. The average times that a pore stays 
open and blocked are 



feonP(l) (feoff + h)9 feoff + fe t ' 

(39) 

respectively. The inverse of the mean time between suc- 
cessive capture events defines the capture rate uj c : 

lo c ps — i— = (feoff + k t )9. (40) 

-tb + -to 

Upon defining the normalized capture rate Q c = cj c /u> ff , 
we formally express 

n c = ( J*) ps ( e -W + ^V) 9, (41) 

where the voltage-dependent expressions for fc D ff and 
fc t have been used and the occupation 9 is de- 
termined as a function of all physical parameters 
(Tfl,/, q, V,w, £z,,Pao,k>oe,a>on, and w t ) by solving Eq. 
34, or using Eqs. 36 or 37. Thus, for fixed applied volt- 
age V, the capture rate as a function of salt ratio E^ is 
proportional to 9, and is plotted in Fig. 5(a). 

Figure 6(a) shows the normalized capture rate Q c as a 
function of voltage V for different salt ratios E^ . For fixed 
El > 1, increasing V initially decreases £l c by virtue of 
the prefactor u> se~f qV . Only at larger V does the term 
L0tV/ijj o e increase the overall capture rate. For small El, 
not only is <£•+ slightly decreased, the (-potential and 
EOF are increased, repelling the analyte away from the 
pore, particularly at high voltages. 

Fig. 6(b) plots the capture rate as a function of voltage 
at various relative translocation rates. For modest Tr 
and El, and relatively small translocation rates wt/w ff, 
the initial increase in £l c with V arises predominantly 
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r„=0.02, aya^O.1. w=0.2 r s =0.1, Z L =2, w=0.2 




voltage bias V voltage bias V 



FIG. 6: (a) Normalized capture rate Q c as a function of bias 
voltage, (a) The effect of varying salt ratio El. (b) The effect 
of varying k on /k t = u on /ujt with fixed El = 2. 

from an increase in small 9 (Fig. 5) which decreases T a 
(Eq. 39) despite the decrease in k Q g. However, for larger 
V, the decrease in k g is not compensated by the slower 
increase in 9 (Fig. 5). Only at large voltages does the 
k t = Lo t V term come into play to increase f2 c linearly. For 
larger translocation rates tu t , the pore is cleared faster by 
annihilation of the analyte into the opposing reservoir, 
preventing the initial increase in 9 with V, as well as 
the initial decrease in T that would increase the overall 
capture rate. 

SUMMARY AND CONCLUSIONS 

We have modeled the underlying electrokinetics to 
quantitatively describe capture of charged analytes by 
nanopores in the presence of salt gradients. Our analytic 
analysis shows that the electrostatic potential near the 
pore mouth, often neglected, can be sufficient to be an im- 
portant determinant in the capture of charged particles. 
We also showed that imposed salt gradients locally change 
solution conductivity locally, altering the potential distri- 
butions. In particular, higher salt (higher conductivity) 
in the non-analyte chamber decreases the voltage drop 
in that reservoir and across the nanopore, increasing the 
analyte blocking probability. As a function of salt ratio, 
our analysis at small analyte concentrations predicts that 
both the blocking probability 9 and the capture rate f2 c 
always increases as Sl increases. The basic mechanism 
provides a physically consistent and testable explanation 
for recently observed increases in capture rate with salt 
ratio T, L > 1 [5]. 

Electroosmotic flow also affects analyte capture. By 
itself, hydrodynamic flow (e.g., from electroosmosis) into 
the analyte reservoir sweeps particles away from the pore, 
dramatically lowering the blocking probability. However, 
when the non-analyte chamber contains a higher salt con- 
centration such that El > 1, the same repelling fluid flow 
can also change the local conductivity structure such that 



the potential $ + felt by the charged analyte at the pore 
mouth initially increases with flow rate. For small pore 
charge/permeability Tr, wc find that repelling EOF ac- 
tually increases the overall attraction of charged analytes, 
particularly when Sl and the effective analyte charge 
q are large. Finally, our analysis shows that the cap- 
ture rate is sensitive to the translocation rate w t . When 
u>t ^> w fi, and nearly all captured particles are anni- 
hilated via translocation into the receiving reservoir, the 
capture rate increases with bias voltage V , except at very 
low V. 

We considered only outer solutions, accurate in regions 
where Q(r) w outside the charged boundary layer at 
the solution-membrane interface. However, since we fo- 
cus on the capture of analyte from the bulk reservoir, the 
main factor is the voltage drop across the pore with other 
electrostatic details within the pore relatively unimpor- 
tant. Although our outer solutions do not hold inside 
pores with small radii a < k^ 1 , the ionic current flow 
across such pores still induces a slowly decaying electro- 
static potential $ r (r) in the bulk analyte chamber, albeit 
with a small amplitude determined by an effective, small 
pore aspect ratio w. The correspondingly smaller po- 
tential $ r (r) in the right chamber would give a smaller 
capture rate at the same analyte density. We expect the 
analyte capture rate by smaller pores to have the same 
functional dependences as in our mathematical model, 
but with a smaller effective w. Note that the EOF veloc- 
ity U is proportional to w 2 (through T(U,T,l) in Eq. 27), 
while $+ scales as w. Therefore, we also expect EOF to 
become less important than direct electrostatic effects for 
sufficiently small pores. 

In addition to the charge-neutral approximation, our 
analysis relies on a number of other assumptions, includ- 
ing an effective analyte charge q and pore surface charge 
a that are independent of the local ionic strength S(r). 
Moreover, we have assumed right-cylindrical pores, that 
the occupation and bulk analyte density can be approxi- 
mated using a mean field assumption, and that the molec- 
ular details of the capture and translocation can be de- 
scribed using simple kinetic rates. Although some of these 
assumptions can be lifted in more detailed models and nu- 
merical analyses (for example, in a model of conical pores 
[28]), our simple model embodies the essential physics of 
the problem and wc expect our results to be qualitatively 
predictive. 
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